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EfRcient methods for the description of the non-Markovian dynamics of open systems play an im- 
portant role in many proposed applications of quantum mechanics. Here we review some of the most 
important tools that are based on the projection operator techniques of nonequilibrium statistical 
mechanics. The standard product-state projection is generalized to a new class of correlated projec- 
tion superoperators that allow the treatment of strong memory effects, and lead to a non-Markovian 
generalization of the Lindblad equation. 



I. INTRODUCTION 



Relaxation and decoherence processes are key features of the dynamics of open quantum systems [l| . In the standard 
^S) ] approach one tries to develop appropriate master equations for the open system's reduced density matrix ps which 
is given by the partial trace taken over the environmental variables coupled to the open system. Invoking the weak- 
coupling assumption one can formulate in many cases of physical interest a Markovian quantum master equation for 
. the reduced density matrix, expressing the dynamical laws for the irreversible motion of the open system. 
I ' However, the theoretical description of quantum mechanical relaxation and decoherence processes often leads to a 
(— I , non-Markovian dynamics which is determined by pronounced memory effects. Strong system-environment couplings 
' 3) 9 ) correlations and entanglement in the initial state (J, Q , interactions with environments at low temperatures 
^ \ and with spin baths @ , finite reservoirs 0, [1] , and transport processes in nano-structures Q can lead to long memory 
^ / times and to a failure of the Markovian approximation. 

' Here, we will review the most important features of a systematic approach to non-Markovian quantum dynamics 
^ which is known as projection operator technique [lol . fill , fl^ . [l3l |. This technique is based on the introduction of a 
^ certain projection superoperator V which acts on the states of the total system. The superoperator V expresses in 
OA , a formal mathematical way the idea of the elimination of degrees of freedom from the complete description of the 
states of the total system. Namely, if p is the full density matrix of the composite system, the projection Vp serves 
to represent a certain approximation of p which leads to a simplified effective description of the dynamics through a 
reduced set of relevant variables. 

With the help of the projection operator techniques one derives closed dynamic equations for the relevant variables 
Vp . We will discuss two different approximation schemes. The first one is based on the Nakajima-Zwanzig equation 
[HI which represents an integrodifferential equation for Vp with a certain memory kernel. The second scheme 
employs a time-convolutionless master equation for Vp, i.e a time-local differential equation with a time-dependent 
generator [3, [l^, [l^, 0, [l^ . These equations are used as starting point for the derivation of effective master 
equations through a systematic perturbation expansion. 
5J] ■ In the standard approach to the dynamics of open systems one chooses a projection superoperator which is defined 
" by the expression Vp — ps® Po, where po is some fixed environmental state. A superoperator of this form projects the 
total state p onto a tensor product state, i. e., onto a state without any statistical correlations between system and 
environment. Many examples for this product-state projection are known in the fields of quantum optics, decoherence, 
quantum Brownian motion, quantum measurement theory, and coherent and optimal quantum control. It is typically 
applicable in the case of weak system-environment couplings. The corresponding perturbation expansion is usually 
restricted to the second order (known as Born approximation), from which one derives, with the help of certain further 
assumptions, a Markovian quantum master equations in Lindblad form (20l. [2ll. [2^ . 

A possible approach to large deviations from Markovian behavior consists in carrying out the perturbation expansion 
to higher orders in the system-environment coupling. However, this approach is often limited by the increasing 
complexity of the resulting equations of motion. Moreover, the perturbation expansion may not converge uniformly in 
time, such that higher orders only improve the quality of the approximation of the short-time behavior, but completely 
fail in the long-time limit [1^. 
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We will discuss here a further strategy for the treatment of highly non-Markovian processes which is based on the 
use of a correlated projection superoperator [2^ [H, [S^, [28l.l2^. By contrast to the product-state projection, a 
correlated projection superoperator projects the total state p onto a system-environment state that contains statistical 
correlations between certain system and environment states. We will discuss a representation theorem for a large class 
of such projections which are appropriate for the application of the projection operator techniques, and develop a 
corresponding non-Markovian generalization of the Lindblad equation. 



II. THE STANDARD PROJECTION OPERATOR METHOD 



We investigate an open quantum system S that is coupled to some environment E. The corresponding Hilbert 
spaces are denoted by Tis and Ti, e , respectively. The state space of the composite system is thus given by the tensor 
product space 

n = Hs®HE- (1) 

The states of the composite system are represented by density matrices p on Ti. satisfying the physical conditions of 
the positivity and the normalization: 

p > 0, trp = 1, (2) 

where tr denotes the trace taken over the total state space TC. The partial traces over Hs and He will be denoted by 
trg and tr^;. 



A. Nakajima-Zwanzig projection operator technique 

A central goal of the theory is to develop efficient strategies for the description of the behavior of the reduced 
density matrix which is determined by the partial trace over the environmental state space, 

PS = t^EP- (3) 

The basic idea of the projection operator techniques is to regard the operation of taking the partial trace over E 
formally as a map V defined by 

Pp= (trsp)®po- (4) 

For a fixed environmental state po this defines a linear transformation which maps any density operator p on the total 
state space H to a density operator Vp on Ti. and has the property of a projection operator: 

r'^ = V. (5) 

Being a map acting on operators, V is often called a projection superoperator. The complementary projection is 
defined by 

Q = I-V, (6) 

I being the identity map. Note that according to Eq. @ the reduced system's state is obtained from the projection 
V p by taking the partial trace over the environment: 

PS = t^EiVp}- (7) 

The Hamiltonian of the composite system is of the form 

H = Ha + Hi, (8) 

where Hq denotes the unperturbed part, usually given by the sum of a system Hamiltonian Hs and an environmental 
Hamiltonian He, and Hj represents the interaction. In many cases it is convenient to formulate the dynamics in the 
interaction picture with respect to Hq in which the density matrix p(t) of the total system is governed by the von 
Neumann equation 

^^p{t)^-^[HJ{t),p{t)]^C{t)p{t). (9) 
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The operator 

Hi{t) = e'""'Hie-'"°' (10) 

represents the Hamiltonian in the interaction picture and C{t) the corresponding Liouville superoperator. 

The Nakajima-Zwanzig (NZ) projection operator technique yields a closed equation of motion for the relevant part 
Vpit) of the density matrix and, hence, for the reduced density matrix ps{t)- To simplify the presentation we assume 
that the condition 

PC{h)C{t2)...C{t2n+l)V^0 (11) 

holds true. This condition is in fact satisfied in many applications. Moreover, we suppose that the initial state 
satisfies Vp{0) = p(0). The projection Vp{t) is then governed by a homogeneous integrodifferential equation, the 
Nakajima-Zwanzig equation: 

j^Vp{t) = j'^dhlC{tM)Vp{h). (12) 

The memory kernel lC{t, ti) is given by 



K.{t,ti) ^VC{t) Texp 



dt2QC{t2) 



QC{ti)P, (13) 



where T denotes the chronological time ordering. 

The memory kernel IC{t, ti) is in general a very complicated superoperator whose determination is in most cases as 
complicated as the solution of the full system's dynamics. Therefore, one usually tries to determine it by a perturbation 
expansion in powers of the strength of the system-environment coupling. The lowest order contribution is given by 
the second-order equation of motion: 

^Vp{t)^ / dtiVC{t)C{ti)Vp{ti). (14) 
dt Jo 

Higher orders are obtained with the help of the general expression (|13p for the memory kernel. 



B. Time-convolutionless projection operator technique 

There exists an alternative expansion technique based on the projection superoperator V which is known as time- 
convolutionless (TCL) projection operator method. By contrast to the NZ approach, the TCL method leads to an 
equation of motion for the relevant part of the density matrix which represents a time-local differential equation of 
the general form 

±Vp{t) = K.{t)'Pp{t). (15) 

Here, K,{t) is a time-dependent superoperator, called the TCL generator. It should be stressed that the TCL equation 
(fT5)l describes non-Markovian dynamics, although it is local in time and does not involve an integration over the 
system's past. In fact, the TCL equation takes into account all memory effects through the explicit time-dependence 
of the generator /C(t). 

To obtain the time-local form of the TCL equation one eliminates the dependence of the future time evolution on 
the system's history through the introduction of the backward propagator into the Nakajima-Zwanzig equation. This 
enables one to express the density matrix at previous times ti < tin terms of the density matrix at time t and to derive 
an exact time-local equation of motion. We remark that the backward propagator and, hence, also the TCL generator 
may not exist, typically at isolated points of the time axis. This may happen for very strong system-environment 
couphngs and/or long integration times; an example is discussed in [ij. 

Again, one can develop a systematic perturbation expansion for the TCL generator which takes the form lC(t) = 
K,2{t)-\-K,4:{t) + . . . The various orders of this expansion can be expressed through the ordered cumulants [30l.l3lll32.l33j 
of the Liouville superoperator C{t). For instance, the contributions of second and fourth order to the TCL generator 
are given by [l| 



lC2{t) = 




dtiVC{t)£{ti)V, 
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and 



dti / dt2 / dt^ 
^0 Jq 



X 



VC{t)C{t2)VCiti)Cih)V - VC{t)C{U)VC{ti)C{t2)V 



In second order the TCL master equation takes the form 



d /■* 

j^Pp{t) = j dtiVC{t)C{h 



(16) 



which should be contrasted to the NZ equation (|T4|) . 

It is important to reahze that the NZ and the TCL technique lead to equations of motion with entirely different 
structures and that, therefore, also the mathematical structure of their solutions are quite different in any given order 
[3^ . It is difficult to formulate general conditions that allow to decide for a given model whether the NZ or the 
TCL approach is more efficient. The assessment of the quality of the approximation obtained generally requires the 
investigation of higher orders of the expansion, or else the comparison with numerical simulations or with certain 
limiting cases that can be treated analytically. It turns out that in many cases the degree of accuracy obtained by 
both methods are of the same order of magnitude. In these cases the TCL approach is of course to be preferred 
because it is technically much simpler to deal with. 

In the NZ equation ([T^ as well as in the TCL equation we made use of the initial condition 'Pp{0) = p(0). 
According to the definition ([4]) of the projection V this condition is equivalent to the assumption that p(0) represents 
an uncorrected tensor product initial state, p(0) = ps{0) ® po- For a correlated initial state one has to add a certain 
inhomogeneity to the right-hand side of the NZ or the TCL equation which involves the initial conditions through the 
complementary projection Qp{0) = (/ — 'P)p(O). A general method for the treatment of such correlated initial states 
within the TCL technique is described in [1]; for a recent study on their influence in weakly coupled systems see also 
Refs. M,M. 



C. Markovian limit and quantum dynamical semigroups 

With the standard projection defined in Eq. the TCL equation is equivalent to the following master 
equation for the reduced density matrix, 

d /■* 

JtPsW = -J^ dtitTE{[Hi{t),[Hi{h),psit)®Po]]}. (17) 

This equation provides an appropriate starting point for an approximation scheme which is known as Born-Markov 
approximation and which eventually leads to a Markovian quantum master equation in Lindblad form 

j^Psit) = ICpsit) 

= -t[Hs,psit)] + Y.(^RxPsml-l{RiRx,Psit)}y (18) 

Here, K. is a time-independent generator, the Lindblad generator, involving a Hermitian operator Hs and arbitrary 
system operators R\. Therefore, it generates state transformations of the form 

^t-Ps{0)^Ps{t), $t-e'=*. (19) 

is called a quantum dynamical map and the set of transformations 

{^t\t>0} 

is referred to as quantum dynamical semigroup. Under certain technical conditions, it can be shown that the form of 
the Lindblad generator guarantees the preservation of the positivity and normalization of the density matrix, as well as 
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the complete positivity of the dynamical transformation Vice versa, any completely positive quantum dynamical 
semigroup has a generator of the form (jlSp . This is the well-known Gorini-Kossakowski-Sudarshan-Lindblad theorem 

The microscopic derivation of the master equation from the TCL equation (fT7|) requires the validity of several 
approximations, the most important one being the so-called Markov approximation. This approximation presupposes 
a rapid decay of the two-point correlation functions of those environmental operators that describe the system- 
environment coupling. More precisely, if te describes the temporal width of these correlations and the relaxation 
time of the system, the Markov approximation demands that 

TE < TR. (20) 

This means that the environmental correlation time te is short compared to the open system's relaxation time tr. 

The Markov approximation is justified in many cases of physical interest. Examples of application are the quantum 
optical master equation describing the interaction of radiation with matter, and the master equation for a test particle 
in a quantum gas [stI. [ssI. [39] . However, strong couplings or interactions with low-temperature reservoirs can lead 
to large correlations resulting in long memory times and in a failure of the Markov approximation. In the following, 
the quantum dynamics of an open system is said to be non-Markovian if the time-evolution of its reduced density 
matrix cannot be described (to the desired degree of accuracy) by means of a closed master equation with a (possibly 
time-dependent) generator in Lindblad form. 

If the two-point environmental correlation functions do not decay rapidly in time the second order of the expansion 
cannot, in general, be expected to give an accurate description of the dynamics. For instance, this situation arises for 
the spin star model discussed in Ref. [1^, where the second-order generator of the master equation increases linearly 
with time such that the Born-Markov approximation simply does not exist. 

More importantly, the standard Markov condition (|20[) alone does not guarantee, in general, that the Markovian 
master equation provides a reasonable description of the dynamics. This situation can occur for finite and/or struc- 
tured reservoirs that cannot be represented by a Bosonic field or a collection of harmonic oscillator modes. In such 
cases a detailed investigation of the influence of higher-order correlations is indispensable in order to judge the quality 
of a given order. The model discussed in Ref. 25] represents an example for which the standard Markov condition 
is satisfied although the expansion based on the projection (j4|) completely fails if one truncates the expansion at any 
finite order. In such cases strong non-Markovian dynamics is induced through the behavior of higher-order correlation 
functions. 

We conclude that in general one can judge the quality of a given projection superoperator and a given expansion 
technique that is based on it only by an investigation of the structure of higher orders. The standard projection and 
the corresponding Lindblad equation are not reliable if higher orders lead to contributions that are not bounded in 
time, signifying the non- uniform convergence of the perturbation expansion [25j |. 

III. CORRELATED PROJECTION SUPEROPERATORS 

The performance of the projection operator techniques depends of course on the properties of the microscopic 
model under study, in particular on the structure of the correlation functions of the model. However, it also depends 
strongly on the choice of the superoperator V. Several extensions of the standard projection ^ and modifications of 
the expansion technique have been proposed in the hterature (see, e.g., Refs. |40l l4ll. |42|). 

The projection defined by Eq. ^ projects any state p onto a tensor product ps (8) po that describes a state without 
statistical correlations between the system and its environment. Here, we introduce a more general class of projection 
superoperators that project onto correlated system-environment states and are therefore able to describe strong 
correlations and non-Markovian effects 1241. 



A. General conditions 

We assume that our new class of maps V represent superoperators with the property of a projection, i.e., = V. 
As a consequence, the whole machinery of the projection operator techniques described in Sec.|TT]can be applied also 
to the new class of correlated maps. 

Within the projection operator techniques the projection Pp should represent a suitable approximation of p. We 
therefore require that for any physical state p the projection Vp is again a physical state, i. e., a positive operator 
with unit trace. This means that "P is a positive and trace preserving map, namely 

p>0=^Vp>0, triVp} = trp. (21) 
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Our class of projection operators is assumed to consist of maps of the following general form, 

V^Is®^. (22) 

Here, Is denotes the unit map acting on operators on Hs, and A is a linear map that transforms operators on He into 
operators on He- A projection superoperator of this form leaves the system S unchanged and acts nontrivially only 
on the variables of the environment E. As a consequence of the positivity of V and of condition (|22p the map A must 
be A^s-positive, where Ns is the dimension of Hs- In the following we use the stronger condition that A is completely 
positive, because completely positive maps allow for a simple mathematical characterization (see Sec. IIIIBI) . 

Let us discuss the physical implications of these conditions. According to Eqs. ([5|) and ([22|) the map A must itself 
be a projection, namely A^ = A. Moreover, since V is trace-preserving, the map A must also be trace-preserving. 
Hence, we find that A represents a completely positive and trace-preserving map (CPT map, or quantum channel) 
which operates on the variables of the environment and has the property of a projection. The action of V may also 
be interpreted as that of a generalized quantum measurement which is carried out on the environment. A further 
physically reasonable consequence of the positivity of A and of Eq. ((22)) is that V maps product states to product states, 
and, more generally, separable (classically correlated) states to separable states. This means that the application of 
■p does not create entanglement between the system and its environment. 

Using Eq. (|22p and the fact that A is trace-preserving we get 

PS = tvEP = trE{Pp}. (23) 

This relation connects the density matrix of the reduced system with the projection of a given state p of the total 
system. It states that, in order to determine ps, we do not really need the full density matrix p, but only its projection 
Pp. Thus, Vp contains the full information needed to reconstruct the reduced system's state. 



B. Representation theorem 



What is the explicit structure of the projection superoperators satisfying the basic conditions formulated above? 
This question is answered by a representation theorem [2J| which states that V fulfills the condition of Sec. IHI Al if 
and only if it can be written in the form 

Pp^J2^''E{A^p}(g,B„ (24) 

i 

where {Ai} and {Bi} are two sets of linear independent Hermitian operators satisfying the relations 

tTE{B,A,} = <5„, (25) 
Y,{trEB^)A, = Ie, (26) 



<B, > 0. (27) 



Equation (|25|) guarantees that is a projection superoperator, Eq. (|26|) ensures that V is trace-preserving, and 
Eq. ([TT]) is equivalent to the condition of complete positivity (T denotes the transposition). 

The standard projection ^ that projects onto an uncorrelated tensor product state is obviously of the form of 
Eq. In fact, if we take a single A — Ie and a single B — po the conditions (PS)) . ([^5]) . and ([77]) are trivially 

satisfied and Eq. ([24]) obviously reduces to Eq. ([4]). Of course, a projection Vp of the form of Eq. (j24p does not in 
general represent a simple product state. We therefore call such V correlated projection superoperators. They project 
onto states that contain statistical correlations between the system S and its environment E. In the following we will 
consider the case that one can find a representation of the projection with positive operators 

A, > 0, B, > 0. (28) 

Equation is then trivially satisfied. Without restriction we may assume that the Bi are normalized to unit trace, 

trEB^ = 1, (29) 

such that condition (pS)) reduces to the simple form 

^A^Ie. (30) 
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Under these conditions V projects any state p onto a state which can be written as a sum of tensor products of positive 
operators. In the theory of entanglement (see, e.g., the recent review [43|]) such states are called separable or classically 
correlated. Using a projection superoperator of this form, one thus tries to approximate the total system's states by 
a classically correlated state. The general representation of Eq. ([M]) includes the case of projection superoperators 
that project onto inseparable, entangled quantum states. We will not pursue here this possibility further, and restrict 
ourselves to positive Ai and Bi in the following. 



C. Correlated initial states 



As mentioned already in Sec. Ill BI a homogeneous NZ or TCL equation of motion presupposes a tensor product 
initial state if one uses the standard projection superoperator However, this is no longer true for the correlated 
projection defined by Eq. ([M]) . In fact, the general condition for the absence of an inhomogeneous term in the NZ 
equation (fT2|) or the TCL equation p5|) is given by 

Pp(0) = p(0). (31) 
According to Eq. this condition is equivalent to the assumption that p(0) takes the form 

p(0) =^p,(0)®B„ (32) 

i 

where 



p,(0) = tr£;{A,p(0)}>0. (33) 

Equation ([5^ represents in general a correlated initial state. Hence, a great advantage of the correlated projection 
superoperators is given by the fact that they allows the treatment of correlated initial states by means of a homogeneous 
NZ or TCL equation j2J]. 



D. Conservation laws 



A crucial step in applications of the correlated projection operator technique is the construction of an appropriate 
projection superoperator V. An important strategy for this construction is to take into account the known conserved 
quantities of the model under study. 

Suppose C is a conserved observable. A good choice for the projection superoperator V will then be a projection 
that leaves invariant the expectation value of C, i.e. that satisfies the relation 

tr{Cp} = t^{C{Vp)}. (34) 

To bring this relation into a more convenient form wc introduce the adjoint V'^ of the projection superoperator V . 
The adjoint map is defined with the help of the Hilbcrt-Schmidt scalar product 

{X,Y)^ix{X^Y] (35) 

for the space of operators acting on the state space of the total system through the relation 

{X,VY) = {V^X,Y). 

This allows us to write Eq. (|34p in the form tr{Cp} = iY{{V'^ C)p^ . Requiring this to hold for all states p we get the 
relation 

V'^C = C. (36) 

The adjoint of the projection p4l) is obtained by interchanging the role of At and Bi. Hence, condition ([36)1 can be 
written explicitly as 

V^C ^Y.^^E{B^C}®A,^C. (37) 

i 

This equation represents a condition for the projection superoperator V on the basis of a known conserved quantity 
of the underlying model. It ensures that the projection superoperator leaves invariant this quantity and that the 
effective description respects the corresponding conservation law. 
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IV. GENERALIZATION OF THE LINDBLAD EQUATION 



Once a projection superoperator has been chosen the projection V p{t) of the time-dependent total system's state 
p{t) is, according to Eq. (|24p . uniquely determined by the dynamical variables 

p^{t)^trE{A,p{t)}. (38) 

To be specific we assume in the following that that the index i takes on the values j = 1, 2, . . . , n. Since we require 
that the Ai are positive, the pi{t) are positive operators. From Eq. (l23l) we find the connection to the reduced density 
matrix, 

Ps{t)=Y,Pr(t), (39) 

i 

and the normalization condition takes the form 

iTsPs{t)^Y.^rsP^{t) = l■ (40) 

i 

Hence, we see that the reduced system's state is uniquely determined by a set of n (unnormalized) density operators 

pi{t). 

Our formulation leads to a natural question, namely what is the analog of the Lindblad equation (|18p in the case 
of a correlated projection superoperator? To answer this question we first observe that the time-evolution leads to a 
transformation of the form 



{p^m {p^m. 



(41) 



transforming any initial set of positive operators Pi(0) > into another set of positive operators Pi{t) > at time 
t > Q. This transformation can conveniently be described with the help of an auxiliary n-dimensional Hilbert space 
C" and a fixed orthonormal basis {|«)} for this space. Then we can identify the collection of densities pi{t) with a 
density matrix Q(t) on the extended space 



Tic 



through the relation 



(42) 



(43) 



This density matrix can be regarded as a block diagonal matrix 



( Pi{t) 

P2{t) 



\ 





Pn{t) j 



(44) 



with blocks pi(t) along the main diagonal. Moreover, the reduced density matrix ps{t") is obtained by the partial trace 
of p(i) taken over the auxiliary space. 

In close analog to Eq. (jl9[) the dynamics may now be viewed as a transformation 



Vi ■■ ^'(0) ^ Q{t) 



(45) 



that preserves the block diagonal structure. It is important to emphasize that Vt is not a quantum dynamical map in 
the usual sense because it is not an operation on the space of states of the reduced system, but rather a map on the 
extended state space. In fact, the transition from q{Q) to the reduced density matrix ps(0) = Pi(0) is connected 
with a loss of information on the initial correlations, such that from the mere knowledge of ps{Q) the dynamical 
behavior cannot be reconstructed. 

It may be shown that Vt can be extended to a completely positive map for operators on Ticxt- Hence, we can 
construct an embedding of the dynamical transformation into a Lindblad dynamics on the extended state space. This 



9 



is achieved by the requirement that there exist a Lindblad generator /C acting on operators of the extended state 
space which preserves the block diagonal structure: 

l^(^P^^\i){^\j =^/C,(pi,...,p„)® (46) 

such that the time-evolution can be represented in the form 

J2 pm = e'^* p,{0) . (47) 

One can show that a Lindblad generator IC with this property exists if and only if the densities pi (t) obey the master 
equation 

^p,it) = -^ [H\p,it)] fep,(t)i?f - i {i?f i?r,ft(t)}) . (48) 

The W are Hermitian operators on Hs, while the R^-^ may be arbitrary operators on Hs- The details of the proof of 
this statement can be found in Ref. [24] . 

Equation ([^5]) represents the desired non-Markovian generalization of the Lindblad equation for the case of 
a classically correlated projection superoperator. This equation has many physical applications. In fact, master 
equations of the form of Eq. have been derived by several authors and applied to various models featuring 
pronounced non-Markovian effects ^ lil HI [13, [H, [H . 



V. CONCLUSIONS 



We have reviewed the theoretical treatment of non-Markovian quantum dynamics within the framework of the 
projection operator techniques. It has been shown that an efficient description of strong non-Markovian effects is 
made possible through the construction of correlated projection superoperators V. The central idea behind this 
construction is to take into account large system-environment correlations by an extension of the set of dynamical 
variables. In fact, employing a correlated projection superoperator instead of a product-state projection, one enlarges 
the set of dynamical variables from the reduced density matrix ps to a collection of densities pi describing system 
states that are correlated with certain environmental states. 

General physical conditions for a large class of correlated projection superoperators have been formulated, demand- 
ing essentially that V can be expressed in terms of a projective quantum channel that operates on the environmental 
variables. These conditions lead to a representation theorem for correlated projection superoperators and to a non- 
Markovian generalization of the Lindblad equation that is capable of modelling long memory times and large initial 
correlations, while preserving the physical conditions of positivity and normalization. 

The method developed here has many applications to physically relevant models featuring non-Markovian dynam- 
ics. The investigated class of projections does not exhaust all possibilities. Future investigations should include 
the formulation of further classes of correlated projections, the study of time-dependent generators, as well as the 
application of correlated maps that project onto nonseparable, entangled quantum states. 
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